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Hydrodynamic simulations are used to make predictions for the integrated elliptic flow coefficient 
i>2 in yfs — 5.5 TeV lead-lead and y/s = 14 TeV proton-proton collisions at the LHC. We predict a 
10% increase in V2 from RHIC to Pb+Pb at LHC, and V2 ~ in p+p collisions unless rj/s < 0.08. 



INTRODUCTION 

Much work has been done recently using viscous hydro- 
dynamics to study the properties of gold-gold collisions 
at the Relativistic Heavy Ion Collider (RHIC) A 
measurement of particular interest is the elliptic flow co- 
efficient V2, the second moment in the azimuthal angle of 
the distribution of emitted particles (cf. @), which allows 
to extract information about material constants (such as 
viscosity) of the high density nuclear matter created at 
RHIC. Using the knowledge gained at RHIC, it should 
be possible to predict experimental results at the Large 
Hadron Collider (LHC), which will collide lead ions at a 
maximum center of mass energy of y/s = 5.5 TeV per nu- 
cleon pair compared to yfs = 200 GeV gold ions at RHIC. 
If experimental data on, e.g., V2 from LHC is close to the 
hydrodynamic model prediction, this would confirm that 
real progress has been made in understanding nuclear 
matter at extreme energy densities; if far away, it may 
indicate that the successful hydrodynamic description of 
experimental data from RHIC was a coincidence. 

Regardless of the outcome, the advent of the RHIC ex- 
periments clearly has lead to major progress in the theory 
and application of hydrodynamics to heavy- ion collisions. 
A few years ago the form of the hydrodynamic equations 
in the presence of shear viscosity rj was still unresolved, 
with different groups keeping some terms while neglect- 
ing others |6fi9|. For the case of approximately confor- 
mal theories, where the viscosity coefficient for bulk — 
but not shear — becomes negligible, all possible terms to 
second order in gradients were derived in Ref. 10| , and 
their relative importance investigated in Ref. [1]. Three 
of the groups performing viscous hydrodynamic simula- 
tions now agree on these terms y, |4j , while another 
group uses a different formalism that gives match- 
ing results. While this development still leaves out the 
consistent treatment of bulk viscosity, the quantitative 
suppression of elliptic flow by shear viscosity is there- 
fore essentially understood. From comparison of viscous 
hydrodynamic simulations to experimental data llj, Il2| , 



one can infer an upper limit of the ratio of shear vis- 
cosity over entropy density, rj/s < 0.5, for the matter 
produced in Au+Au collisions at yfs = 200 GeV [l|, 
comprising extractions by other methods 13- HI . A size- 
able uncertainty for this limit comes from the fact that 



the initial conditions for the hydrodynamic evolution are 
poorly known, with the two main models, the Glauber 
and Color-Glass-Condensate (CGC) models, giving dif- 
ferent results for the elliptic flow coefficient [1]. This 
difference can be understood to originate from the dif- 
ferent initial spatial eccentricity e x in the Glauber/CGC 
models. The eccentricity is defined as 



(y 2 -x 2 ) 

(y 2 + x 2 ) 



(1) 



where the symbols () denote averaging over the initial 
energy density in the transverse plane, e(x,y). 

Indeed, it had been suggested [16| that the elliptic flow 
coefficient vi at the end of the hydrodynamic evolution 
would be strictly proportional to the initial spatial eccen- 
tricity, V2/&X °c const., if the fluid was evolving without 
any viscous stresses for an infinitely long time. This is to 
be contrasted with experimental data indicating a pro- 
portionality factor of total multiplicity over overlap area 
V2/&X oc dN/dY/S 'overlap ,16]. Total multiplicity here 
refers to the total number of observed particles N per 
unit rapidity V, while the overlap area is calculated as 



^overlap = 71 V ( x2 ) (j/ 2 > • 



(2) 



Ideal fluid dynamics does not adequately describe the 
latest stage of a heavy-ion collision (the hadron gas) , be- 
cause of the large viscosity coefficient in this stage [l7j . 
Therefore, the hydrodynamic stage lasts only for a fi- 
nite time (e.g. until all fluid cells have cooled below the 
decoupling temperature), resulting in a dependence of 
v-ije-x on dN/dY. Also, viscous effects affect the propor- 
tionality between Vi and e x , leading to a behavior that 
is qualitatively similar to that observed in the data [J]. 

One of the objectives of this work is to extend the en- 
ergy range for fluid dynamic results of V2 / e x from Au+Au 
collisions at top RHIC to Pb+Pb collisions at top LHC 
energies, as well as to study the dependence on shear vis- 
cosity. If in the future either e x or the mean rj/s becomes 
known, these results can thus be used to constrain the 
respective other quantity from experimental data. On 
the other hand, the values of shear viscosity for which 
the Glauber/CGC models match to experimental data 
at top RHIC energies have been extracted in Ref. 0, HH 
for Au+Au collisions. Since rj/s averaged over the system 
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evolution is not expected to be dramatically different for 
Pb+Pb collisions at the LHC, another objective of this 
work is to obtain a prediction for the elliptic flow coeffi- 
cient for the LHC based on the best- fit values to RHIC. 

Finally, the feasibility of detecting elliptic flow in p+p 
collisions at s/s = 14 TeV at the LHC is being dis- 
cussed [l9| . As a reference for other approaches and ex- 
periment, it interesting to study the possible size and 
viscosity dependence of v-i under the hypothetical as- 
sumption that the bulk evolution following p+p collisions 
could be captured by fluid dynamics. 



Beam 


Initial cond. 


. dY , 


Ti [GeV] [GeV] 


r [fm/c] 


Gold 


Glauber 


800 


0.34 


200 


1 


Gold 


CGC 


800 


0.31 


200 


1 


Lead 


Glauber 


1800 


0.42 


5500 


1 


Lead 


CGC 


1800 


0.39 


5500 


1 


Protons 


Glauber 


6 


0.400 


14000 


0.5 


Protons 


Glauber 


6 


0.305 


14000 


1 


Protons 


Glauber 


6 


0.270 


14000 


2 



TABLE I: Central collision parameters used for the viscous 
hydrodynamics simulations (T/ = 0.14 GeV for all). 



SETUP 

To make predictions for nuclear collisions at LHC en- 
ergies, we use our hydrodynamic model that successfully 
described experimental data at RHIC [H, HH and make 
modifications to the input parameters appropriate for the 
higher collision energies at the LHC. 

As a reminder, the hydrodynamic model 1] is based 
on the conservation of the energy momentum tensor (loj 



IP 1 " = rjV^u 
Ai 



m 



A, 



2r? 2 



2r/ 



where e,p and u M are the energy density, pressure, 
and fluid 4-velocity, respectively. D = u^D^ and 
V Q = A^D M are time-like and space-like projec- 
tions of the covariant derivative D^, where A' Il/ = 
q^f _ anc [ we remind the compact notations 

A {ll B v) = (A^+A^-fA^A^)^ and 
cj m „ = | (Vj/Uju — VhUjj). For relativistic nuclear col- 
lisions it is convenient to follow Bjorken [2(| and use 
Milne coordinates proper time r = \/t 2 — z 2 and space- 
time rapidity £ = atanh-|, in which the metric becomes 
9V = diag(l, — 1, — 1, — t 2 ), and assume that close to 
£ = 0, the hydrodynamic degrees of freedom are approx- 
imately boost-invariant (£ ~ Y). 

The hydrodynamic equations D^T^ = then consti- 
tute an initial value problem in proper time and trans- 
verse space, and are solved numerically (see Ref. [H). The 
input parameters for hydrodynamic evolution are the 
equation of state p = p(e) and the first (second) order hy- 
drodynamic transport coefficients rj (m, Ai, A2, A3). The 
values for Xi } 2,3 have been found to hardly affect the 
boost-invariant hydrodynamic evolution for Au+Au col- 
lisions at RHIC [l|, so here they are generally set to zero. 

The mechanisms leading to thermalization (the onset 
of hydrodynamic behavior) are not well understood in 
nuclear collisions. Therefore, it is not known how the 
thermalization time tq at which hydrodynamic evolution 
is started depends on the collision energy. Barring further 



insight, we start hydrodynamic evolution for the LHC at 
the same time as for the RHIC simulations (to = 1 fm/c). 
At this time, the initial conditions for the transverse en- 
ergy density e(x,y) are given by the Glauber or CGC 
model, respectively, the fluid velocities are assumed to 
vanish, and the shear tensor IP"' is set to zero (other 
values for IP 1 " do not seem to affect the final results 
[HHH). F° r brevity, we refer to Ref. [H for the details of 
the Glauber and CGC models, but for the expert reader 
note that we use the Woods-Saxon parameters of radius 
Rq = 6.4 (6.6) fm and skin depth \ — 0-54 (0.55) fm for 
gold (lead), and assume a nucleon-nucleon cross section 
of a = 40 (60) mb for ^ = 200 (5500) GeV collisions. 

The overall normalization of the initial energy density 
(parametrized by the initial temperature at the center Tj) 
was adjusted to match the experimentally observed mul- 
tiplicity at RHIC; by analogy, for LHC the normalization 
is adjusted to match the predicted multiplicity [22J-125J. 
Since we lack detailed knowledge about its temperature 
dependence, the ratio of shear viscosity to entropy den- 
sity rj/s is set to be constant during the hydrodynamic 
evolution (equal to the average over the spacetime evolu- 
tion of the system) . The relaxation time coefficient rn is 
expected [HI HH to lie in the range ^(e + p) ~ 2.6 — 6. 
The equation of state (EoS) can in principle be provided 
by lattice QCD. While at present there are points of dis- 
agreement between lattice groups about, e.g., the precise 
location of the QCD phase transition, there is consen- 

Therefore, 



sus that it is an analytic crossover [27, 
we use a lattice-inspired EoS [29( that is consistent with 
both the current consensus and perturbative QCD; also, 
since it resembles 27J, we expect that using a different 



lattice EoS will have a minor effect on our results. 

Once a given fluid cell has cooled down to the decou- 
pling temperature 7/, its energy and momentum are con- 
verted into particle degrees of freedom using the Cooper- 
Frye freeze-out prescription [30]. A value of Tf = 0.14 
GeV was determined by matching to RHIC data and will 
also be used for LHC energies, assuming that it is mostly 
determined by local conditions, and less so by initial en- 
ergy density, system size or collision energy. The distri- 
bution of the particle degrees of freedom may be further 
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FIG. 1: (Color online) Anisotropy ([3]) divided by (JTJ, as a 
function of initial entropy (Q divided by ((2j. Shown are 
results from hydrodynamic simulations for y/s = 200 GeV 
Au+Au (RHIC) and y/s = 5.5 TeV Pb+Pb collisions (LHC). 
For comparison, experimental data for V2 from RHIC [3l[, di- 
vided by e x from two models [l3|], is shown as a function of 



measured 



divided by ([2]). See text for details. 



evolved using a hadronic cascade code (as in Ref. [33|). 
or in a more simple approach the unstable particle res- 
onances are allowed to decay, without further evolving 
the stable particle distributions. In both cases, the total 
multiplicity and particle correlations (such as the elliptic 
flow coefficient) are then calculated from the stable par- 
ticle distribution (cf. [1]). Surprisingly, it was found in 
Ref. [l], Q that the momentum integrated elliptic flow co- 
efficient for charged hadrons — to good approximation — is 
equal to half the momentum anisotropy, 



V2 



1 

2 6p 



1 J dxdy t xx - ryy 

2 J dxdy T xx + Tvv 



(3) 



Since the momentum anisotropy is a property of the fluid, 
it is independent on the details of the freeze-out proce- 
dure and only mildly dependent on the choices of To, T/. 
Unlike at RHIC where pairs of To and Tf could be fine- 
tuned to fit the particle spectra at central collisions, no 
such extra information is available for the LHC. Hence 
Eq. ([3]) may provide the most reliable way of determining 
the elliptic flow of charged hadrons, and will be used in 
the following. Similarly, one can use the total entropy 
per unit spacetime rapidity ^| in the fluid as a proxy for 
the total (charged hadron) multiplicity per unit rapidity 



W (^y 1 ) with a proportionality factor (HHi] 



dS_ 



dS_ 

dY 



4.87 



_dN 
dY 



7.85- 



dN, 



ch 



dY 



(4) 



Note that for a gas of massive hadrons in thermal equi- 
librium at Tf — 0.14 GeV the ratio of entropy to par- 
ticle density is ~ 6.41, but the decay of unstable reso- 
nances produces additional entropy, resulting in Eq. 



FIG. 2: (Color online) Anisotropy (|3]) prediction for y/s = 5.5 
TeV Pb+Pb collisions (LHC), as a function of centrality. Pre- 
diction is based on values of rj/s for the Glauber/CGC model 
that matched y/s = 200 GeV Au+Au collision data from 
PHOBOS at RHIC ([Hj], shown for comparison). The shaded 
band corresponds to the estimated uncertainty in our pre- 
diction from additional systematic effects: using e p /2 rather 
than^ (5%) Q; using a lattice EoS from J2EJ] rather tha n [27 
(5%); not including hadronic cascade afterburner (5%) [3S 



Since results from RHIC suggest there is only approxi- 
mately 10% viscous entropy production during the hy- 
drodynamic phase 0, HiJ j the entropy 4p at r = To can 
be used to estimate the final particle multiplicity. In 
the case of the LHC, the world average for the predicted 
charged hadron multiplicity for central Pb+Pb collisions 
at y/s = 5.5 TeV [13, ^f- ^ 1800, can be used to esti- 
mate the total entropy at r = r , and hence the overall 



normalization Tj of the initial energy density (see Tab.[T]). 

Using Eqs. (|3l4p for the multiplicity and elliptic flow 
allows to make predictions for the LHC without having 
to model the hadronic freeze-out, which should make the 
results more robust. However, as a consequence one does 
not get information about the momentum dependence of 
the elliptic flow coefficient, prohibitin g de tailed compar- 
ison with predictions by other groups 23l |37|. 



RESULTS 

With the initial energy density distribution fixed at to, 
the hydrodynamic model then gives predictions for the 
ratio of V2/e x at the LHC. In Fig.Q] the results are shown 
for three different values of shear viscosity, for two dif- 
ferent initial conditions and two different beams/collision 
energies (Au+Au at y/s = 200 GeV, Pb+Pb at y/s = 5.5 
TeV). The resulting values for V2/e x seem to be quasi- 
universal functions of the total multiplicity scaled by the 
overlap area Overlap, only depending on the value of rj/s 
(and, to a lesser extent, the collision energy). The devi- 
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ations of the RHIC simulations from the universal curve 
can be argued to arise from a combination of the fi- 
nite lifetime of the hydrodynamic phase at s/s = 200 
GeV and the presence of the QCD phase transition, and 
is strongest for ideal hydrodynamics, in agreement with 
earlier findings 

Also shown in Fig. [T] is experimental data for the el- 
liptic flow coefficient for Au+Au collisions at RHIC, nor- 
malized by e x from a Monte-Carlo calculation (inclu ding 
fluctuations) in Glauber and CGC models (see Ref. [13[ 
for details). Since e x is not directly measurable, the dif- 
ferently normalized data gives an estimate of the overall 
size of Vije x at RHIC. Directly matching experimental 
data on v 2 using a hydrodynamic model with an initial 
e x specified by the Glauber or CGC model, a reason- 
able fit was achieved for a mean value of rj/s ~ 0.08 and 
■q/s ~ 0.16, respectively [H, [Hj]. Under the assumption 
that the average 77/s is similar for collisions at RHIC and 
the LHC (along with the assumptions discussed above), 
one can make a prediction for the integrated elliptic flow 
coefficient for charged hadrons as a function of impact 
parameter (or more customarily the number of partici- 
pants Ap ar t, cf. [l|). The result is shown in Fig. [2] As 
can be seen, we expect integrated v 2 at the LHC to be 
about ten percent larger than at RHIC, which is less than 
the prediction by ideal hydrodynamics [1^] , and inagree- 
ment with the extrapolations by Drescher et al. [23| . 

Finally, using the charge density parametrization of the 
proton p(b) in Ref. [Ic| as an equivalent of the nuclear 
thickness function in the Glauber model (cf. [l|) one ob- 
tains an estimate for the shape of the transverse energy 
density following a relativistic p+p collision. Using the 



predicted multiplicity at mid-rapidity ^ 6 [24j, |25j for 



y/s = 14 TeV p+p collisions at the LHC, one can again 
use Eq. (j4|) to infer the overall normalization of the energy 
density (or Ti) at r = tq (see Tab.H]). As a "Gedankenex- 
pcriment" one can then ask how much elliptic flow would 
be generated in LHC p+p collisions if the subsequent evo- 
lution was well approximated by boost-invariant viscous 
hydrodynamics. One finds that for ideal hydrodynamics 
Ie. V2 0.035 for integrated \v 2 \ in minimum bias col- 
lisions (cf. (28) in [l[), while for rj/s > 0.08, v 2 typically 
changes by almost 100 percent when varying the relax- 
ation time ^(e + p) between 2.6 and 6 and varying To by 
a factor of two. This indicates that for rj/s > 0.08, the 
hydrodynamic gradient expansion does not converge and 
as a consequence it is unlikely that elliptic flow develops 
in p+p collisions at top LHC energies. If experiments 
find a non- vanishing value for integrated 1 1^2 1 > 0.02 in 
minimum bias p+p collisions, this would be an indication 
for an extremely small viscosity -q/s < 0.08 in deconfined 
nuclear matter. 

To conclude, viscous hydrodynamics can be used to 
make predictions for the ratio of v 2 /e x as a function 
of multiplicity and rj/s. Assuming a multiplicity of 
"' x " ~ 1800 for the matter created in Pb+Pb collisions 



at LHC, as well as rj/s similar to RHIC, we predict the 
integrated elliptic flow for charged hadrons to be 10% 
larger at the LHC than at RHIC. We expect v 2 measure- 
ments in p+p collisions to be consistent with zero, unless 
the shear viscosity is extremely small {rj/s < 0.08). 
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